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I relate some coefficients encountered when computing the functional 
determinants on spheres to the central differentials of nothing. In 
doing this I use some historic works, in particular, transcribing the el- 
egant symbolic formalism of Jeffery (1861) into central difference form 
which has computational advantages for Euler numbers, as discovered 
by Shovelton (1915). I derive sum rules for these, and for the central 
differentials, the proof of which involves an interesting expression for 
powers of sech x as multiple derivatives. I present a more general, 
symbolic treatment of central difference calculus which allows known, 
and unknown, things to be obtained in an elegant and compact fashion 
gaining, at no cost, the expansion of the powers of the inverse sinh, 
a basic central function. Systematic use is made of the operator 2 
asinh(D/2). Umbral calculus is employed to compress the operator 
formalism. For example the orthogonality/completeness of the facto- 
rial numbers, of the first and second kinds, translates, umbrally, to 
T(t(x))=x. The classic expansions of multiple angle cosh and sinh in 
terms of powers of sinh are thence obtained with minimal effort. 
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1. Introduction 

In a work, [1] , involving quantum field theory calculations of the effective action 
on spheres (as the logarithm of the operator determinant of the field propagation 
operator), in order to show the equivalence of two explicit expressions for this 
'logdet', I was obliged to look at an expansion of the powers of the sech function. 
Some details are in the following section. 

As background, the power series expansion of secx goes back to Euler, the 
coefficients being named Euler numbers in the 1820's to 1840's, e.g. by Scherk and 
Raabe. They were, and still are, called secant numbers. Their evaluation is distinct 
from those of tan x which are related to the Bernoulli numbers. The series expansion 
of powers of sech a; yields generalised Euler numbers. 

The study of Bernoulli numbers, their relatives and their generalisations re- 
mains strong after nearly three hundred years. For this reason, it is difficult to 
know whether any displayed result is really new and, not unexpectedly, the expan- 
sion in [1] turns out to be known and I wish here to go into more detail about this 
interesting topic which has involved me in an analysis of the existing, and sometimes 
ancient, literature. The Dilcher online bibliography contains most (not all) of these 
references. Whether anyone other than the historically minded has actually used, 
or even examined, this old material is open to question. 

The relevant calculus turns out to be that of central differences and in section 5 
I transcribe some of the symbolic formalism given by Jeffery, [2] , who uses forward 
differences, into central calculus. This allows various sum rules and identities to be 
derived, which might be new, as well as some (known) series expansions. 

In sections 6 and 7 I further develope the symbolic approach to central calculus 
in a more general way presenting some standard and non-standard applications in 
a (possibly) novel fashion. In section 8, I combine this approach with the umbral 
notation, generating functions, and other things, emerging effortlessly. 

2. The expansion and central factorials 

The integral needed in [1] is, 




(1) 



and, in order to allow partial integration, one requires the expansion, 




p=0 
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for odd powers in terms of even derivatives. 2 The coefficients satisfy the recursion, 
(difference equation), 

g; = (2r - if g;- 1 + g;zl , (3) 

together with Q r _ x = = an d the initial value Qq = 1. Special values are 

g r r = 1. They are directly related to central difference factorial numbers (e.g. 
Steffensen, [3]) as is easily seen by repeated action of the basic, and old, relation, 

sech n x = —. ^ (D 2 -in- 2) 2 ) sech n - 2 x , D = ^- , 

(n - l)(n -2) K v ' ' ' dx 

which yields, 

sech 2r+1 x = { j^r{D 2 - (2r - f) 2 ) {D 2 - (2r - 3) 2 ) . . . (D 2 - 1 ) sechx 
(2r). 



[ 2r +l]-l 

2 2r (— J sechx, 



(4) 



in Steffensen's notation. The coefficients of the expansion of the central factorial 
in ordinary powers of D are, by definition, the central factorial numbers of the 
first kind (or, up to a factor, the 'central differentials of nothing') 3 . The complete 
iteration, (4), shows that these are equal to the sum of the squares of the first r 
odd integers taken a certain number of times (c/Ely [4], Jeffery [2]) but I will not 
continue with this combinatorial view just now. 

The precise connection between the coefficients in (2) and the central differen- 
tials of nothing is, 4 

nr _ ( -.\r-pp2r ( 2r ) ! D 2p+1 n [2r+l] 

P ~ { } (2p+l)! (5) 

= (-1 ) r - p 2 2r (2r)! t(2r + 1, 2p + 1) , 

where I have introduced the notation of Riordan, [6] , for the central factorial num- 
bers, t. The recursion, (3), is then equivalent to those in Steffensen, [3], and in [6] 
and [7]. 

Writing (2) this way round was the most convenient for the purposes in [1]. 
One could, of course, invert it and, in this form, a related expansion is given by 



2 I have improved the definition of the coefficients in order to make them integers. 

3 One might call them 'central Stirling numbers of the first kind'. 

4 A slightly more complicated treatment, but essentially the same, is given by Zhang, [5]. 
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Stern, [8], which I wish to relate to (2). He actually gives two expansions, one for 
even derivatives, 

^sechx = (-l) p sechx^(-l) fc 4 2p) tanh 2fc a:, (6) 

and one for odd, 
d 2 P +1 



1 sechx = (-l) p+1 cosechx V(-l) fc a[ 2p+1) tanh 2fc+2 x 
dx 2 P +1 ^ k 

k=0 

p 

= (-l) p+1 sechx tanhx ^(-l) fc a^ 2p+1) tanh 2fc x . 



(7) 



fc=0 



One sees the appearance of 'derivative polynomials' which have been subject to 
more recent investigations e.g. [9]. 

Stern gives the coupled recursions, 

a^ +1) = (2*+l)a< 2 »> + (2fc + 2)a< 2 ;> 
a^ +2) =2ka^ 1) + (2k + l)a^ +1) . 

The coefficients Q r can be connected to the by substituting (6) into (2) 
giving, 

r p 

(1 - tanh 2 x) r = J2 g p 5Z(- 1 ) fca i 2p) tanh2fc x > 

p=0 fc=0 

whence, 

j^G;a^ = {2r)\Q. (9) 

p=s 

Similarly, substituting (7) into the derivative of (2) produces 

j^G;af^ = {2r + l)\Q. (10) 

p=s 

The second recurrence in (8) is equivalent, directly, to the binomial recurrence 

r \ r — s /r\ /r\ 
s + 1 ) ~ s + 1 \s) ' V0/ _ 
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3. Generalised Euler numbers 



Generalised Euler numbers can be denned by the power series (I use Norlund's 
notation, [10,11]), 

The original Euler numbers correspond to n = 1, E 2v = E 2 J . 

Setting x to zero in (6) and using the MacLaurin expansion, yields 

— a — a > 

where the second equality follows on differentiating (7), and so Stern refers to the 
other coefficients, a£, as Euler numbers of higher order but this name is more often 
used for the numbers defined by (11). Stern's term reappears in some recent works, 
e.g. [9]. 

By differentiating (2) 2u times, one can easily obtain an expression for the odd 
generalised numbers in terms of the ordinary ones, 

4l r+1) = j^XyrG r p E 2u+2pJ Vv. (12) 

This relation can be found in Ely, [4], where the coefficients, Q, are expressed in the 
combinatorial fashion. 

Schwatt, [12], pp. 68-71, has an extensive algebraic treatment of expansions of 
powers of sechx and also of the derivatives of these in the Stern form, (6), (7). 

One could look upon (12) as a means of computing the (odd) generalised num- 
bers since the recursion (3) for Q is easily solved and the first order Euler numbers 
are fairly well studied. In particular an efficient recursion method is adopted by 
Knuth and Buckholtz, [13], which it is worth recalling and then comparing with 
earlier work. (See also Brent and Harvey, [14]. ) 



4. Computing Euler numbers 



Knuth and Buckholtz use trigonometric functions and, by differentiation, de- 
duce the relation, 

d n 

secz = sec z(E nQ + E nl tanz + E n2 tan 2 z + . . .) , (13) 



dz 
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and the recursion, 



E n +i,k = kE njk _ x + (k + l)E njk+1 , E ok = 5 0k . (14) 

Since E n = E n0 this allows an efficient computation of the Euler numbers. 5 

The earlier Stern recursions are the same as (13) which is easily seen by splitting 
(13) into two by the parity behaviour under z — > —z. Then one sees that, 

E 2p,2k = a l P > E 2p,2k+1 = > 

and 

F _ _(2p+l) F _ n 

^2p+l,2k+l — a k ' ^2p+l,2k — u - 

Stern uses the recursions to derive many arithmetical properties, as do Knuth 
and Buckholtz, who were unaware of Stern's work. 

Historical computations of the secant numbers were performed by Euler and 
by Scherk (in 1825). I do not give these continental references as they can be 
found in the 1893 text by Saalschiitz, [15], as well as a list of the first 14 numbers. 
Glaisher extended the list to 27 and Joffe to 50. The effort seems to be more of a 
computational challenge than necessity. 

Herschel, [16,17], (1816 and 1820) gave finite difference expressions suitable for 
numerical computation (in terms of forward differences of nothing) and also related 
Euler numbers to Bernoulli numbers via tangent numbers. 

Herschel's theorem applied to sechx yields, 

^''-''• (hahW 1 ^ 1 (15) 

which is manipulated in various ways. However, a more elegant, and practical, 
formulation was found by Shovelton, [18], in terms of the central difference operator, 
5, as, 6 



1 r\2n — I i\n/\ 1 X2 , 1 r4 , ,;2n \ n 2n 



E 2n = (" 1 ) n rT ^ = (- 1 ) 71 - 2 5 + ¥ 5 + • • • + V" 6 ' ° ' (16) 



^ A simple algorithm implementing (13) is given by Brent and Harvey. Programming this in 
DERIVE (not a fast CAS) I calculated the first 1000 Euler numbers in 24 seconds using an 
Athlon 610e processor. The values could be used in a look-up table for fast access. 

6 Shovelton developes his results via the Bernoulli polynomials (the 'Bernoullian function') and 
introduces central differences as a simplifying tactic. 
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which can be obtained directly from (15) by noting that 5 = A E 1 / 2 (E being the 
translation operator, E = 1 + A) so that (e.g. [18], [3] pp. 10, 11, [19]), 



S = E 1 ! 2 — E~ x l 2 , i.e. 2 + 5 2 = E + E' 1 , 



and also, 



(£1/2 + £-1/2)2 =4 + 5 2_ 



(17) 



After use of the equivalence f(E n ) s = n s f(E) s , Vn (Herschel, [20] VII 
§(8).) equations (17) and (15) allow the alternative compact expression, also given 
by Shovelton, [18], 



(-1)' 



1 



^2n 22n ( 1 + £2/4)1/2 

-h\8 







In 



_ (-1)' 



2 2r 



1 + 



1 2 2 



^7 + 



24+- + 



-^5 



v2n 



2 2r 







In 



(18) 



(See Butzer et a/, [7].) 

Joffe, [21], undertook the numerical evaluation of E 2n using (16) in terms of 
the central differences of nothing which were computed using the recursion formula. 
An early table of values is provided by Shovelton, [18]. 

The expansion of the powers of sec is formally no harder and by the same 
process one has the coefficients, [18], 



E. 



(r) 
2n 



(-1)' 



1 



(l + 5 2 /2Y 
1 ( —r 







2n 



1 + 



A 2 ■ 1 
1 J 6 + ¥ 



5 4 + ...+ 



n 



r2n 







In 



(19) 



which has the advantage over (12) of displaying the polynomial nature of E^J in r. 
The form corresponding to (18) is, 



E {r) - 



(-iy 



2 2n (l + ( y2/4)r/2 







2n 



2 2n 



1 + 



5 1 

2 4 + 



+ 



n 



S 2n 
2^ 



(20) 







2n 



The central differences of nothing are, to a factor, the central factorial numbers 
of the second kind, e.g. Steffensen, [3], and [6] and [7]. 
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Earlier expansions of sech n x exist but are generally expressed in terms of the 
forward differences of nothing. For example Jeffery, [22] , gives the neat looking, but 
computationally somewhat awkward, form, 



V 2 + A / ' ' 

An expression, essentially equivalent to (19), has been given more recently by 
Liu and Zhang, [23]. 

5. Jeffery's theorem, central differences and Euler relations 

The early works of Herschel, Jeffery, Boole and others are characterised by the 
extensive use of symbolic methods. Shovelton does likewise, but in terms, partly, 
of central differences. The 1909 textbook by Thiele, [24], is a rare example of their 
application by a continental worker. Another is the classic text by Steffensen, [3], 
both geared to interpolation questions. There is also an 1865 work by Hansen 
that I have not been able to access. Bickley, [19], gives a useful summary, without 
references, of standard operator relations. Chapter 6 in Riordan, [6], is a good, 
modern pedagogic treatment. 

The elegant papers by Jeffery contain applications of four basic theorems con- 
cerning the transformation of symbolic expressions. (Herschel's much earlier work 
also involves a lot of these manipulations.) The formalism is expressed in terms of 
forward differences. Shovelton's results show that considerable simplification ensues 
if central differences are brought into play and I wish here to transcribe some of Jef- 
fery's formulae into central form and apply them to obtain Euler number relations 
as an illustration. 

For this purpose, I return to the alternative formula, (18), for the Euler num- 
bers, and write it equivalently as, 



Rather than expanding the operator in powers of 5 2 directly (as in (18)), I 
substitute for the 2n in central factorials using, 




(-l) n 2 







2n 



2 2n £1/2 + E -l/2 
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to give the intermediate formula, 



2 2n ^ 2 (2v) 



(21) 



I am following Jeffery, [22] §3, and, as there, a theorem is employed to re- 
express the action of the pre-operator, here, sech.D/2. The theorem is the central 
difference form of Theorem D in Jeffery, [2] , and reads (I prove this in section 6) , 

/(£>) N = HM. Qln-l] = . . . = < EUSP1 1 (22 ) 

v ' d5 d5 n v ' 

which can be run left to right or right to left. Applied to powers of sechD/2, taking 
it all the way to the right trivially gives the expansion (20), and so nothing new. 

More interesting is to work right to left. For example, set n = 2v + 1 and 
choose, 



f D f D D 

f(D) = / dS sech — = / dD cosh — sech 



2 (23) 



= D, 

(since 5 = 2sinh D/2) so that (21) becomes, 



/ i\n n z2vr\2n 

^ = ^E O0l2 " +1] -7^r. (^) 

which is identical to (18) as the value of the particular central differential in (24) 
is standard, [3] (see below). This result is given in [7] but only by comparing the 
actual values. 

Another conclusion can be drawn from the identity, 

sec h^0^]= J D0^+ 1 ], (25) 
by comparing with the expansion (11) to give the sum rule for Euler numbers, 



V E2n D 2n 0^ = D 0[ 2s+1 l = (~^\ 
^ (2n)!2 2 - 2 2s 1 s J 



1.3 (2s - 1) 



2 



(26) 



which might be novel. Given the central differentials, it could be used to compute 
the Euler numbers, but there is no need. 
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Related manipulations can be performed for the generalised Euler numbers, 
(20), although the simple relation, (25), becomes more complicated. As an example, 
in order to further illustrate the theorem, (22), I look at the numbers. Instead 
of (23), f(D) in (22) is given by, 

f(D) = I dD cosh — sech 3 — , 
J J 2 2 (27) 

= 2tanhL>/2, 

which results in, 

sec h 3 ^0^ = 2tanh^0^ +1 ]. (28) 

The tanh could be expanded in order to evaluate the right-hand side but it is better 
to apply (22) for a second time by choosing, 

D , D 

— tanh — 

2 2 (29) 

= 4coshL>/2, 
yielding the successive equivalences, 



dD cosh ^ tanh ^ 



sech 3 - 0^1 =2tanh-0^ +1 ] = 4 cosh - 0^+ 2 l , (30) 
2 2 2 v 1 

the last of which can now easily be evaluated as a finite sum of central differentials, 

n "+1 n 2fc n [2^+2] 

4 cosh - 0^ +2 ] = 4 V - 

4cosn 2 u 4^ 2 2k { 2k)\ ' 

generalising (25). 

Spelling out the details, the resulting expression for the generalised Euler num- 
bers follows as (21), 



^2n " 22n 2^ SeCh 2 " {2v)\ 

(31) 



n ~ 2 ^ ^ 2 2fc (2/c)! (2i/)! ' 



2 2 



which is considerably more complicated than the direct expansion, (20), but allows 
us, therefrom, to infer the (even) sum rules, 

^2 fcQ[ 2, +2] ( 2 „)!/-f\ 



^ 2 2fc -2(2A;)! 2 2 ^ 
9 



which can be confirmed numerically. 

As another example, I give the (odd) sum rule that results from looking at 

p(5) 

4 ^ D 2 fe+ l [2, + 3] _ (2v)\(-\\ 

3^ (2fc + l)! - 2*» U J' (33J 
the formula corresponding to (30) being, after three integrations, 

sec h 5 -0^ = -sinh J D0^+ 3 ]. (34) 
2 3 

The general case depends on integrating the relation, 7 

-P^f-TT^-V^™. r = l,2,..., (35) 



r (2r — 1)!! \(isinhx 

which I cannot find listed. I need not write out the corresponding sum rules as 
examples (32) and (33) are sufficient. 

It is clear that there are a large number of sum rules derived by using the 
different transformations, such as (30) and (34), arising from (22) and by comput- 
ing higher Euler numbers, E^ +1 \ They would seem to have only a curiosity or 
checking value. I have not seen them elsewhere. 

With (35) I have essentially returned to the direct expansion, (20), since the 
multiple-angle cosh and sinh can be expanded in powers of sinh x in a familiar way, 
the coefficients usually being derived by recursion or by schoolbook trigonometry 
(e.g. Loney, [25], Bromwich, [26], Workman, [27], Schwatt, [12] p.124). 

The expansions can actually be combined neatly to the unfamiliar form, 

e rx = J22 } sinh 5 a: , reZ, (36) 

where I have written the coefficients in terms of central factorials, for future use in 
section 8. (See [7], §3.1.) 

Eq.(35), after some brief algebra, coincides with the direct expansion, (cf 
(20)), 8 

sech 2r+1 a; = >J I 2 I sinh 2s x , sinh a; < 1 , 

s =o V s / 

cf also Schwatt, [12] p. 70. 9 Questions of convergence have been mostly ignored, 
and, in any case, for the finite series involved in (20) are irrelevant. 



7 Constants of integration give zero. 

8 This could be considered a proof of (35), given the expansions, (36), or these could be derived, 
by continuation. 

9 Schwatt also refers to Shovelton. 
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6. Central difference theorems 



(37) 



I develope some basic concepts in central difference calculus using a mainly 
symbolic approach. I do this by first proving the theorem (22) and then passing on 
to other notions and some expansion formulae. Of the few discussions of this topic, 
that by Riordan, [6] , is most valuable. The more recent review by Butzer et al, [7] , 
is also useful. 

The proof of (22) is not obvious so I expose it here, following the discussion by 
Jeffery, [2], for forward difference calculus. For comparison I write out the theorems 
as given by him, 

A : /(A)O n =n^^O n - 1 = ... = n!^^l 

B: /(A)0" =<«*>»-> = ...= *^*>1 
JK ' dD dD n 

C : f(D) 0<»> = (J*"" 1 ) = . . . = n!^l 

D : /(D)0W = ^)0M = ... = «1. 
JK ' dA dA n 

I will not prove them but only say that A and C express the fact that the 
symbols A and D equally commute when acting upon x at x = as they do when 
x is current. 10 

In proving D, Jeffery makes essential use of the operator £ = log(l + D) which 
is a sort of dual, or functional inverse to the relation D = e A — 1 in the sense that, 

e log(l+D)_ 1 = D _ 

I start from the very basic property of differences as the action on factorials, 

F(S) [r] = F(D) r , (38) 
expressed in terms of central quantities for which a basic relation is, 

5 = 2sinh — , i.e. D = 2sinh" 1 -. (39) 
I can therefore define a new function / by the functional relation, 
F(S) = F(2 sinh D/2) = f(D) = f{2 sinh" 1 5/2) , 



10 As noted by Jeffery, A and B appear as Misc. Ex. (1) on p. 242 of Boole, [28]. Since this was 
a Senate House Examination problem, it seems likely it was due to Herschel. 
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and just replace F(5) by f(D) on the left-hand side of (38) while, on the right-hand 
side, replace 5 by D in /, regarded as a function of 5. So, corresponding to if I 
define the operator by 

£ = 2sinh- 1 ^. (40) 
I obtain the handy reciprocal result, 11 



/(D)0H = /(£) 0' 



(41) 



I have not seen this approach in the literature. 

The proof of the central version of Theorem D now follows more or less directly. 
Applying Theorem B, which holds generally, on appropriate interpretation of /(A), 
one finds, 



/(D)0W = /(€) r 

= no i + 

= /'(C) 



(1 + Z> 2 /4)V2 







r-l 



-h\D 



r_1 + 



1 J 2 2 

l 

2 

1 



+ 



-±\ D 4 
24 



+ 







r-l 



1 \ Q r-3 



2 2 



+ 



- \0 



\r— 5 



2 4 



+ 



(42) 



/'(£) 







[r-l] 



/'(£>) 



(l + 5 2 /4)V : 



+ 







2 » u 



r-l] 



2 2 
df(D) 
dd 



+ 



I\ n[ r - 5 l 

2 » u 



2 4 



+ 



Q[r-l] . 



This is equation (22). 

Continuing with some central formalism, I display the expansion of the central 
factorial, 



x 



x 



„£"0 r 



„=i ^ (43) 

e x? n , 

using (41) to rewrite the coefficients and also formally summing the series (1.0 n = 
for n^O). (Compare Jeffery, [2], §5.) 

Furthermore, from the explicit definition of (40), the expansion of the inverse 
sinh function is contained in, 



(2sinh- 1 |) n = re 



n x 



(44) 



11 Jeffery derives his formula in a different way, but I find the present route easier. 
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showing rapidly that the coefficients are just the central factorial numbers of the 
first kind (cf [7] §4.2 for a more expansive treatment.) 

The same expansion is derived by Butzer et al, [7] §4.1, but less attractively in 
my view as the proof involves manipulating several series. See also Charalambides, 
[29] and the classic Schwatt, [12], pp. 6, 7, 125, who credits Cauchy. 

As a further illustration of the symbolic formalism, I derive an intrinsic relation 
between the first and second kind central factorial numbers, i.e. between the £ m n 
and S r s . 

I start from the definition of the second kind, 



5"0 n A „*„..8 V Q 



n 



V\ * — ' v\ 



and equate powers of x to get the compact statement of 'orthogonality', 

fe 0{ O n =n\S^ . (45) 

The 5 on the right-hand side is a Kronecker delta. Equation (45) is really a state- 
ment about inverses, and tantamount to completeness (of polynomial bases), the 
reciprocal relation being, 

8 n e°to m = m\8™ . (46) 

cf Riordan, [6] p.213. 

A more formal expression of completeness is, 



e°- s e°-t = 1 ( 
e°^e os = h 



(47) 



L S 7 

with, 

/(*) h = /(C) , f(8) h = f(8) ■ 
Two further useful relations are 

f(D)O r = f(£)e°- 5 O r (48) 

and 

m = f{D)e -*. (49) 

The first is a corollary of (45) and the second a trivial consequence of the basic 
derivative, D m n = ml S^. 
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Corresponding to (45) and (46) there are two important expansions connecting 
multiple derivatives and multiple differences. I derive these for central differences 
following the forward difference treatment of Boole, [28] p. 24, which is an application 
of Maclaurin's theorem in its secondary form, 

0(f) = <f>{D) e o t , (50) 

where D acts on 0. From the identities, (39), 5 is a function of D and so, setting 
t = D and 4> = 5 n , 

5 n = 5 n e°- D , (51) 

where the S on the right-hand side acts on 0. This is the required expression in 
symbolic form. Expansion of the exponential gives a sum of powers of D and, acting 
on a function of x, yields the concrete representation, 12 

6 n u(x) = — T Dm < x ) • ( 52 ) 

m=n 

Inversely, consider <p = D m as a function of 6, and so set t = 5. According to 
(50), 5 in D then has to be replaced by D. This gives the £ operator, (40), and so 
one has the required symbolic expression, 

D m = re oi , (53) 
yielding the concrete, Newton series form, 

D m u(x) = ^f" S m <*) > 

n=m 

which could be derived from (52) using (45). These expansions are standard, see, 
e.g., [7] §6.1. 

I note in passing that the explicit forms of the operators, (39), (40), are not 
required when deriving the symbolic equations. Writing the relation, (39), as 5 = 
^(D), £ is defined by the functional inverse, £ = ^~ 1 (D) = ^~ 2 (d). A choice of 
*(*) other than (40), or log(l + *), would entail a factorial different from the usual 
ones. 

12 I ignore questions of convergence and assume that the function, /, possesses the qualities 
required for all the equations to exist. 
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7. Expansions 



As Jeffery points out, knowledge of the factorial numbers is useful in many 
expansions. For example, 

/ ^2 sinh -1 (sinh -1 = /^sinh" 1 e x 

= /^2sinh- 1 ^ e°^e 0x 



(54) 



/ 2 

/(do^ + o^ 



where (49) and (43) have been used to give the second and fourth lines. This process 
can be continued. 

As a calculated example, consider the simplest case of (54) which is 



2 sinh" 1 ( sum" 1 -) = £ 0^ • - t + £ 0^ ~ . . . . (55) 
V 2 J 1! 2! 

and require the coefficient of x r ' /r\, 

s=0 S ' 

(This is really going back to the third line of (54). A simple iteration of (44) also 
yields the same result.) The summation index s must be odd, and so therefore must 
r (as one knows). Hence the coefficient of x 2k+1 /(2k + 1)! is 

k £n 2 i+ 1 

p n[2fc+l] _ 5 U . +2j + l 2fe+l 

^ ^(2J + 1)! 5 

- V" (~^\ 1 t2j + l n 2fc+l 

which can be evaluated straightforwardly. 
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8. Using Representative Notation 

A simpler case to (54) follows either as an extension of (44) or directly, 

/ (2 sinh" 1 |) = m e° * = £ 1^- x r , (56) 

r=0 

and allows me to bring in Blissard's seductive 'Representative Notation' (or classical 
umbral calculus, [30]) by writing the primary form of Maclaurin's theorem as, 

where the power f n is to be replaced by the derivative, /^(O), after the expansion 
of the exponential. Then the elegant result, 

/ (2sinh-'|)=^i 1 f^ 

r=Q 

_ ~ /w r 
- 2^ r \ x ' 

r=0 

follows on using (43). One could write this result more enigmatically as, 

f(t)e°-* = ef*, (57) 

where now the power, f r , first represents the central factorial, f^ r \ in which the 
powers of / are still to be replaced by the derivatives, /( n )(0), as above. 
In the form, 

oo f r r i 

/M = £^(2sinh§)'\ 

this expansion is derived, in a complicated way, by Teixeira, [31], using existing 
trigonometric expansions, cf (36). He sets y = i6 and splits the formula into real 
and imaginary parts. The derivation here, in contrast, is rapid and from finite 
difference fundamentals, which is how the trigonometric expansions can be, and 
are, more logically obtained, [32]. A more elaborate, conventional discussion can be 
found in [7] §5. 

Riordan employs the umbral 'method' and I use his notation, (see Eq.(5)), for 
the factorial numbers, which has operational advantages). My object is to combine 
umbral techniques with the Jeffery operator formalism (in central form). 
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I have already introduced the representative, t, by, 

n 

t n (x) = t n (x) = *(™> v) x v = x [n] = e x? n , (58) 
so that, by definition, 

x n = e tS Q n^ f = f ( x } _ (59) 

I now define its reciprocal, T, by, cf [6], 

n 

T n {x) = T n (x) = T(n, u) x» = e xS n . (60) 

On the right I have given the Jeffery central equivalents which enable me to obtain, 
by summation, the umbral expressions, 

e yt = e xt e o y ( 61 ) 

e yT = e x5 e Oy ^ (62) 

where t = t(x) and T = T(x). 

From the definitions of the operators £ and 5, the generating functions, (61) 
and (62), take the explicit forms, cf (56), 

e Vt(x) _ e 2xsinh _1 y/2 ^ j-gg^ 

and 

e yT(x) _ e 2xsinhy/2 ^ 

both of which Riordan (see also Carlitz and Riordan, [33]) derives by series combi- 
natorial manipulation but which here follow almost painlessly. 

The two expressions are related by reciprocity as can be shown in the following 

way. 

I first derive the umbral statement of reciprocity. Umbrally, from Eq.(62) 
by setting x — > t(x) one arrives at the formally neat expression of orthogonal- 
ity / completeness , 

T(£(*)) = *. (65) 

(This also easily results from (60) followed by (58), which gives T n (t(x)) = x n .) 
The reciprocal relation follows on iteration of (65), 

f(T(f(*))) =*(*), (66) 
17 



or, invoking completeness, 

t(T(*))=*. (67) 

One might therefore formally set T = ete. Eqs. (65) and (67) are the umbral 
equivalents of the explicit (45) and (46) . 

Now, making the replacement x — > T(x) and employing (67) turns (63) into 
(64) as promised. 

I emphasise the rather interesting fact that Eq.(63) is identical with the classical 
trigonometric expansion, (36), which was expressed in a relevant way for this reason. 
Adjusting notation, (63), under y — > 2sinhx and x — » r/2, becomes, neatly, 

e rx = e 2tsinhx , t = t(r/2), (68) 

which is the umbral representation of (36) . A trigonometric expansion has therefore 
been obtained with very little effort. 13 

Finally I make a formal functional generalisation in accordance with the re- 
marks at the end of section 6. If £ = \I/ _1 (-D), the generating functions are 

e yt(x) =e x*- 1 (y) j (gg) 

and 

e yT( x ) = e ^(y) ? (70) 
with t and T now standing for mutually inverse generic 'factorial' numbers. 



9. Conclusion 



The compactness and rapidity of the central difference symbolic calculus has, 
I think, been demonstrated. I hope also to have shown that there is a considerable 
body of ignored historical work that could usefully be exploited. Whether the 
particular results are interesting is, of course, a matter of personal taste. Certainly, 
the symbolic approach has an intrinsic elegance and power. I draw attention to the 
convenience of the operator, As shown in the last section it allows the double 
generating functions of Carlitz and Riordan to be obtained with ease and thence 
some classic trigonometric expansions, likewise. 

13 Eq.(36) is to be found as Eq.(28) on p. 215 of [6] derived in a combinatorial way. (Riordan refers 
to Polya and Szego.) The approach could also be considered as an alternative derivation of the 
expansions (36). Riordan also exhibits (68). (To make the connection, set x — > z/2, r — > 2x and 
t(r/2) -> t(x).) 
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